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The ATLAS Collaboration recently reported strong modifications of dijet properties in heavy ion 
collisions. In this work, we discuss to what extent these first data constrain already the microscopic 
mechanism underlying jet quenching. Simple kinematic arguments lead us to identify a frequency 
collimation mechanism via which the medium efficiently trims away the soft components of the jet 
parton shower. Through this mechanism, the observed dijet asymmetry can be accomodated with 
values oi qL that lie in the expected order of magnitude. 



1. Introduction. 'Jet quenching' refers most gener- 
ally to a picture in which high-p^ partons produced in 
hard processes lose, in addition to their vacuum fragmen- 
tation, a significant amount of energy in interactions with 
the QCD matter produced in heavy ion collisions. Mod- 
els implementing such medium-induced parton energy 
loss can account for characteristic features of the data 
measured at the Relativistic Heavy Ion Collider (RHIC) 
over the last decade. In particular, the assumption of a 
significant medium-modification of the final state parton 
shower is supported by the strong suppression of all high- 
Pt hadron spectra by a factor ~ 5 in central heavy ion 
collisions compared to p-p collisions; by the absence of a 
comparably strong medium effect in the photon spectra; 
and by the dependence of the suppression pattern on the 
centrality of the collision. For a more complete account 
of RHIC data, and for an overview of open issues in the- 
ory and modeling of jet quenching, we refer to several 
recent reviews [T-^. 

Almost all the experimental support for jet quenching 
at RHIC comes from the study of leading single inclusive 
hadron spectra and high-p^ triggered leading hadron cor- 
relation functions. However, to test and further constrain 
the basic view of jet quenching as a medium- modified 
final state parton shower, it is not sufficient to charac- 
terize the energy degradation of the leading fragment in 
the shower. One needs also to characterize the pattern 
of energy redistribution amongst fragments without ex- 
perimental biases on the fragmentation pattern. That is 
to say by identifying jet fragmentation patterns in heavy 
ion collisions not via trigger particles but via calorimet- 
ric jet measurements. Motivated by these considerations, 
first calorimetric jet measurements at RHIC have become 
available in the last two years [THS]- However, the max- 
imal center-of-mass energy ^snn = 200 GeV in Au-Au 
collisions at RHIC has limited, so far, these studies to 
jet energies Et < 40 GeV. Moreover, the small cross sec- 
tions for such jets at RHIC clearly limit the possibility 
to perform detailed studies. 

Only a few weeks after the start of the LHC heavy 
ion programme, the ATLAS collaboration released first 
data on the medium-induced modification of calorimet- 



rically measured jets in Pb-Pb collisions at ^/s = 2.76 
TeV [10^ (see also preliminary data from the CMS col- 
laboration [H]). In these heavy ion collisions at the TeV 
scale, one finds already in the first 1.7 iih~^ of data abun- 
dant samples of jets with Et > 100 GeV, supporting the 
expectation that jet samples with Et > 300 GeV are fi- 
nally in experimental reach for runs at design luminosity. 
Most importantly, the first ATLAS data indicate that for 
jets of (9(100) GeV, the modifications of jet fragmenta- 
tion in heavy ion collisions are much more significant than 
the experimental uncertainties that may arise from iden- 
tifying jets in a high- multiplicity environment. This is 
the very basis for a qualitatively novel access to under- 
standing the microscopic dynamics of parton propagation 
in dense QCD matter. 

In the present rapid communication, we present a qual- 
itative discussion of how the first ATLAS data already 
can contribute to constrain our dynamical picture of jet 
quenching. We start with a short discussion of the data. 

2. Discussion of the ATLAS measurements on 
jet quenching. ATLAS analyzed event samples from 
Pb-Pb collisions at ^/sNN = 2.76 TeV in which the trans- 
verse energy of the most energetic jet (leading jet) was 
Et-^ > 100 GeV and where a recoiling jet with transverse 
energy Et2 > ET^min = 25 GeV could be found at az- 
imuthal separation A<p = |^i — ^2! > 7t/2. Both jet en- 
ergies are identified with the anti-/cT algorithm [12^ with 
radius R = 0.4. Since the shape of jets reconstructed by 
this algorithm tends to be close to circular in the r]-(j)- 
plane, we shall loosely refer to the region assigned to a 
jet as a cone of radius R. 

In Fig. [1] (top) we have replotted the measured dis- 
tribution of such events for the 10 % most central heavy 
ion collisions as a function of the fraction x of the leading 
jet energy carried by the recoiling jet, x = Et^/Et^. In 
comparison to p-p collisions, ATLAS observes in Pb-Pb 
collisions a strong reduction of the number of recoiling 
jets carrying a large fraction {x > 0.6, say) of the maxi- 
mally available jet energy within a cone of radius R = 0.4, 
see Fig. [TJtop). In contrast, the number of recoiling jets 
carrying a smaller energy fraction (x < 0.5, say) is en- 
hanced in Pb-Pb collisions. ATLAS also shows data on 



2 



the angular dijet distribution in A0, see Fig. [TJbottom). 
This distribution displays a significant broadening in the 
range it/2 < Acj) < Sir/A that contains a very small frac- 
tion (9(5%) of ah dijets. In the region 37r/4 < A(/) < tt, 
where more than 95 % of all dijets lie, the distribution 
dN/dA(p broadens at best rather mildly. 

Jet finding in the high multiplicity environment of a 
heavy ion collision poses challenges beyond those present 
in p-p collisions. This is the case since in the 10% most 
central Pb-Pb collisions there are of the order of 1500 
[13] charged particles with typical transverse momentum 
larger than 500 MeV which result in (9(100 GeV) back- 
ground energy in any cone of radius R = 0.4. While 
this background energy can be subtracted, background 
fiuctuations could lead to large artificial asymmetries in 
the reconstructed dijet energy. However, the ATLAS 
and CMS collaborations have put such concerns largely 
to rest by embedding simulated dijets without medium- 
modifications both into a Monte Carlo model of the back- 
ground [10 and into real event samples [TT]. This pro- 
vides strong support to the crucial assumption that the 
measured asymmetries are due to medium modifications 
of the jet properties. 

The data shown by ATLAS are normalized to the 
'number of events' and not to the total number A^tot of 
leading jets with energy Et-^ > 100 GeV. For the study of 
jet quenching, one would prefer that yields are normal- 
ized to A/tot- To the extent to which there are leading jets, 
for which no recoiling jet with > Et, min is found by 
the jet finding algorithm, Aevt will be smaller than A^tot- 
Any medium-induced softening of the leading jet to val- 
ues < ET^min (but also any multiplicity-dependent 
inefficiency in finding jet recoils of low Et2 > ET,min 
GeV) will result in a reduction of A^evt that is more pro- 
nounced in Pb-Pb than in p-p. 

These and other considerations prompt us to limit the 
following discussion largely to qualitative statements. We 
start by shortly recalling the qualitative picture on which 
models of jet quenching are based, before interpreting the 
ATLAS data in the context of this picture. 

3. A qualitative picture of dijet production in 
heavy ion collisions. Fig. [2] illustrates how a dijet is 
embedded in a heavy ion collision. In the plane trans- 
verse to the beam direction, a central Pb-Pb collision 
fills a region of diameter 10 fm with dense QCD mat- 
ter. In contrast, the primary partonic process of a dijet 
event in this region occurs at a large momentum transfer 
of 0[Et) and is therefore localized on a point-like scale 
^ 1/Et within the QCD matter. This sharp localiza- 
tion, illustrated by the red dot in Fig. [2j implies that 
typical soft momentum components of the surrounding 
QCD matter cannot resolve the primary hard partonic 
interaction and therefore will not modify it. However, 
the partons produced in the primary hard process may 
traverse a significant path length within the QCD mat- 
ter, and it is during this final state propagation that the 
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FIG. 1: Top: Dijet distributions measured by ATLAS as de- 
scribed in the text, and plotted as a function of x = Et^/Et-y 
for Etx > 100 GeV. Data from proton-proton collisions at 
y/s = 7 TeV are compared to data from the 10 % most central 
Pb+Pb collisions at ^snn = 2.76 TeV. Bottom: The same 
dijet sample, plotted as a function of the azimuthal angle Acj) 
between leading and recoiling jet. 



medium can modify the jet structure. 

We note that already in proton-proton collisions, there 
are characteristic differences between the leading jet and 
its recoil. In particular, requiring a maximal jet energy 
Et-y within a cone of R = 0.4, one selects jet fragmen- 
tation patterns that deposit more than the average jet 
energy fraction inside the subcone of size R = 0.4. In 
the presence of medium effects, further trigger biases can 
occur. In particular, if there is a medium- induced mech- 
anism that degrades the jet energy fraction in a sub-cone 
as a function of in-medium path length then the leading 
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FIG. 2: Schematic illustration of the spatial embedding of a 
dijet event in a heavy ion collision. In a central Pb+Pb col- 
lision, the overlap of the lead ions in the plane transverse to 
the beam direction fills a region of more than 10 fm diame- 
ter with dense QCD matter. The leading jet and its recoil 
propagate through this matter on the way to the detectors. 
Compared to typical time scales in QCD, (9(10 fm/c) is a 
very long time scale for interactions between a set of partonic 
projectiles and the surrounding QCD matter. This allows for 
strong medium-modifications of jets in heavy ion collisions. 

jet will be oriented preferably in a direction in which its 
path length is minimal. This results in a surface bias 
of the location of dijet events entering the experimen- 
tal data sample. On average, the recoiling jet will see 
a larger in-medium path length and will hence suffer a 
more significant medium modification than the leading 
jet. On the other hand, there may also exist a signif- 
icant contribution of dijet events produced tangentially 
to the nuclear overlap region, for which the in-medium 
path length of the recoiling jet is comparable to that 
of the leading jet {corona effect). A quantitative under- 
standing of the medium modification of dijet events must 
thus rely on a detailed geometric modeling of the nuclear 
overlap region and its strong dynamical evolution. 

In the context of this general picture, we emphasize 
some qualitative features in the data: That Fig. [l] (top) 
indicates a significant medium modification of recoiling 
jets supports the argument that there is a significant sur- 



face bias effect. Therefore, we expect that the structure 
of the triggered jets in Pb-Pb collision is comparable to 
that of triggered jets in p-p collisions (this expectation 
is, of course, testable experimentally). The data of Fig. [l] 
(top) imply then that a significant fraction of the recoil 
energy typically contained in a cone of radius R = OA 
has been transported out of it. And Fig. [l] (bottom) im- 
plies that this transverse redistribution of jet energy has 
been accompanied, for most of the dijets, by a relatively 
mild change in the azimuthal correlation of the center of 
the recoiling jet. Therefore, any dynamical modeling of 
the ATLAS data must account for a significant medium- 
induced transverse broadening of the multiplicity distri- 
bution inside the recoiling parton shower as to displace 
an additional fraction of the jet energy away from the 
recoiling jet cone. But this broadening of the multiplic- 
ity distribution must be accompanied by a rather mild 
medium-induced transverse broadening of the distribu- 
tion of recoiling jet centers. 

We now turn to estimates of how much energy of the 
recoiling jet has been carried out of the jet cone, and to 
what extent the first ATLAS data constrain the dynam- 
ical mechanism underlying this redistribution. 

4. Some lessons from data. We are interested in 
estimating by how much the fraction of the energy radi- 
ated outside the cone of the recoiling jet increases from 
p-p to central Pb-Pb collisions. We start by calculating 
the mean value x of the dijet distributions in Fig. [l] (top), 

X = — — [ dxx ^— . (1) 

Nevt J dx ^ ^ 

In general, the mean x will differ from the average frac- 
tional energy (x) = AE{R)/Et^ radiated outside the 
cone of the recoiling jets. Differences between x and {x) 
arise both because the distribution of Fig. [T] is not nor- 
malized by the number of leading jets, and it does not 
include recoiling jets with energy Ej^^ <C -£/T,min' 

Both 

these effects imply that (x) < x. Therefore, we use the 
calculated mean as a lower bound 

{x)pp<0.Q7, {x)p,p,<OM. (2) 

The difference between x and (x) discussed above is more 
pronounced in Pb-Pb than in p-p collisions. Therefore, 
the discrepancy between {x)^^ and Xpp should be smaller 
than that between the values extracted for central Pb-Pb 
collisions, and {x)^^ — {x)p^p^ > 0.13. 

Here, the average fractional energy {x) was normalized 
to the energy Et-^ of the leading jet, reconstructed with a 
fixed cone R = OA. This is generally smaller than the to- 
tal jet energy, Et > Et^. However, as already explained, 
leading jets tend to maximize the energy within the cone, 
which makes Et-^ a good proxy for Et- From jet shape 
measurements [14] we estimate that for Et-^ > 100 GeV 
we have 1 > Et^/Et > 0.8, and 

^>0.8((x)^,-(x)p,p,) >0.1. (3) 
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So, for total energies Et in excess of 100 GeV, medium 
effects in Pb-Pb collisions deposit at least 10 GeV ad- 
ditional energy outside the cone R=0.4 of the recoiling 
jet. 

The estimate Eq. ([3| is only a lower bound on the en- 
ergy lost by out-of-cone radiation in the medium since 
we are assuming that all jets interact with the medium. 
But as discussed above, there is the corona effect: due 
to geometry, a significant fraction of dijet events are pro- 
duced tangentially to the collision region and will differ 
only mildly from dijets in p-p. To give an upper bound 
on (A£^) /Et^ we take account of the corona effect by 
an all-or-nothing mechanism, according to which only a 
fraction {1 — a) of recoil jets interacts with the medium 
while the rest remain unchanged. From the ratios of the 
yields at x ^ 1 in Fig. [l] we estimate that the fraction of 
unmodified jets a is a ^ 0.5. This implies a mean loss of 
energy in the medium 
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< 
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< 0.2. 



a 



(4) 



We conclude from the above that the interaction of the 
jet with the medium must be able to transport at least 20 
GeV outside the jet cone. Therefore, while the suppres- 
sion of single particle spectra at the LHC [16] requires 
a strong degradation of the longitudinal distribution of 
jet fragments, the ATLAS data imply a strong medium 
induced transverse broadening. 

A mechanism that leads both to large energy degra- 
dation and significant broadening is the medium induced 
radiation of a single hard parton at a large angle. How- 
ever, since a significant fraction of the jets lose more than 
20% of their energy, such a mechanism would imply an 
azimuthal displacement of those jets by Acj) > 0.1 lead- 
ing to a clearly distinguishable modification of the az- 
imuthal distribution with respect to p-p. This is in clear 
contradiction to the data in Fig.[l] Thus, the constraints 
imposed by these data are indicative of a mechanism of 
energy loss involving the out of cone emission of many 
soft partons. 

5. The medium as a frequency collimator. The 

picture that energy loss is carried by soft components is 
also supported by simple kinematic arguments. In the 
dense QCD matter produced in heavy ion collisions, par- 
tonic jet fragments interact within a path length of sev- 
eral fm with several constituents of the medium. We 
expect that these interactions show the main qualitative 
features of partonic scattering cross sections, namely they 
are dominated by small-angle scattering and they have a 
relatively weak dependence on center of mass energy. As 
a consequence of the dominance of multiple small angle 
scattering, all components formed in the parton shower 
undergo Brownian motion. Let us denote by q the aver- 
age squared transverse momentum that a partonic com- 
ponent accumulates per unit path-length by this Brown- 
ian motion. In a medium of length L, all jet fragments 



will accumulate on average a transverse momentum ^/q^. 
Therefore, the softest jet fragments of energy 

CO < ^/4L (5) 

will be decorrelated completely from the initial jet di- 
rection. In other words, the medium acts as an efficient 
frequency collimator which trims away the soft compo- 
nents of the jet. 

The picture advocated above and sketched in Fig.|2]is a 
classical one: all partons in the parton shower accumulate 
on average the same transverse momentum due to mul- 
tiple scattering in the medium. The less energetic these 
components are, the more the medium-induced scatter- 
ing will decorrelate them from the jet axis. To make this 
picture quantum mechanically consistent, one requires 
an argument that the softest partonic components in the 
parton shower can be regarded as independent quanta 
while propagating through the medium; that means, soft 
components must decohere from the parton shower on a 
sufficiently short time scale. This point can be argued 
on the basis of several QCD-based calculations of parton 
energy loss that support the following qualitative conclu- 
sion: The typical formation time r, needed for a parton 
to decohere and become an independent quanta, is set 
by the inverse of its transverse energy r ~ p- (measured 

with respect to the direction of its parent parton). In 
the dense QCD medium, the squared transverse momen- 
tum of the gluon grows by Brownian motion, (/c^) ~ ^t, 
and the average formation time for partons of energy uj 
becomes 



(r) 



(6) 



So, in the QCD medium, softer gluons decohere from 
the parton shower at earlier time scales than harder ones 
because medium-induced scattering decorrelates them 
faster from the direction of the parent parton. 

The kinematical argument outlined above shows that 
soft jet components are easily transported out of a cone of 
finite radius. We now turn to the question whether these 
soft components can carry a sufficiently large fraction of 
the total jet energy outside the jet cone. Our discussion 
will be based solely on the mean energy distribution of 
partonic components inside a jet; fluctuations around this 
mean remain to be considered in a more detailed study. 

For a parton shower evolving in the vacuum, intra-jet 
multiplicity distributions as a function of the longitudi- 
nal momentum fraction z ~ uj/Et can be calculated for 
sufficiently small values z within the modified leading log- 
arithmic approximation (MLLA). These distributions are 
functions of both the jet energy scale Et and the resolu- 
tion scale Qf down to which components in the parton 
shower are counted. The single inclusive parton distri- 
bution as a function of ^ = log 1/z obtained by MLLA 
evolution from an initial hard scale Qo = Et sm{R) down 
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FIG. 3: Dashed: Inclusive parton distribution in a jet of en- 
ergy Et — 100 GeV obtained by MLLA evolution from an 
initial scale Qo = Et sin(i?), {R = 0.4) up to a final partonic 
scale Qf — 1 GeV. Solid: same distribution obtained by a 
medium modified MLLA kernel ^5]. Note that due to kine- 
matical constraints, there are no gluons with z < Et sinO/Q. 



to a final scale Qf is shown in Fig. [s] (dashed line). Con- 
sistent with the praxis in Monte Carlo event generator we 
choose Qf = 1 GeV as the lowest scale for partonic evolu- 
tion. At this scale, one sees that high energy jets contain 
already many soft partons. From the single inclusive dis- 
tribution, the average energy fraction of the jet carried 
by partons with energy fraction smaller than z is given 

by 



E{z) 



log 1/z 



:dD 



(7) 



which is shown in Fig. |4] (dashed line). As long as the 
integration involves only soft components, z <C 1, MLLA 
provides a good approximation for E{z)/Et- 

If the frequency collimation of soft partons is the sole 
medium modification, we can estimate qL hy determin- 
ing from Fig. [4] the value z for which the mean fractional 
energy coincides with the bounds on energy loss, Eq. (|3| 



and Eq. (|4|). Since partons with energy o;^ 
are lost from the cone, we obtain 



! Tp2 



< qL 



35 



Ej" 



< qL < 85 



Ej" 



GeV^ (estimate 1) , 



(8) 

with Et the jet energy given in units of = 100 GeV. 

A generic feature of all jet quenching models is the 
enhancement of small z fragmentation partons as a con- 
sequence of medium induced gluon radiation. This effect 
is in particular necessary for describing the strong sup- 
pression of single inclusive hadron spectra measured at 
RHIC and the LHC Ll6j. Bound Eq. ^ neglects this 
effect, since it is based on a vacuum fragmentation func- 
tion. By supplementing the MLLA framework with a 
medium-induced enhancement of parton branching, one 
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FIG. 4: Fraction of the total jet energy, Et, carried by par- 
tons of energies less than cj = zEt obtained from Fig. [3] via 
Eq. ^ in vacuum (dashed) and both in medium (solid) . The 
plot does not extend to z ^ 1 since MLLA is only valid at 
small z. 



can obtain simple models for the longitudinal softening 
of jet fragmentation functions [15]. An example of such a 
medium-enhancement which is roughly consistent with a 
factor 5 suppression of leading hadron spectra, is shown 
in Fig. [3] and Fig. [4] (solid lines). Since there is a larger 
fraction of the total jet energy stored in soft components, 
a collimation up to the same frequency ^/qL leads to a 
larger energy loss. In this case, we obtain 



30 



Eo) 



< qL < 60 



Ej" 



GeV^ (estimate 2) . 



(9) 

These estimates are subject to various uncertainties. 
Amongst the model-intrinsic ones, we mention the choice 
of the final resolution scale Qf that has significant im- 
pact on the distributions shown in Fig.jsjand Fig. |4] (solid 
lines). Moreover, there are harder partonic components 
of order O (^^/qL/li) that will be partially moved out- 
side the jet cone. Including these components properly 
will require a discussion of fluctuations. Taking the above 
estimates at face value, an extraction of q demands infor- 
mation about the distribution of in- medium path length. 
To arrive at first reference values, we note that L ~ 6 
fm (L - 10 fm) yields 5 < ^ < 10 GeVVfm (3 < g < 
6 GeV^/fm). It is clear, that these first estimates can- 
not replace detailed model studies that must account for 
both the suppression of single inclusive hadron spectra 
[16] and the quenching of true jets [lOl [TT] . 

6. Conclusions. In the discussion of RHIC data on 
single inclusive hadron suppression, emphasis was placed 
on the strong longitudinal softening of the hardest parton 
in the shower. For this, radiative parton energy loss was 
identified as the dominant mechanism. However, there 
was little experimental constraint so far on the angular 
distribution of this radiation. As a consequence, models 
of radiative parton energy loss pT-^ did not focus on an 
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accurate treatement of transverse broadening. Neverthe- 
less, some exploratory studies of jet broadening have been 
performed in the framework of radiative energy loss [17]- 
[19]. 

Here, we have stressed that trimming soft components 
away from the jet parton shower via frequency collima- 
tion provides an efficient mechanism for reducing the en- 
ergy in the jet cone. The main dynamical components 
invoked in this picture have been implemented in various 
models of jet quenching. In particular, collisional parton 
energy loss was identified in some model studies as a sig- 
nificant subleading mechanism [20]. Also, the analysis of 
medium responses to jet propagation led to the proposal 
of mechanisms that contribute to the transverse broad- 
ening of jet fragmentation in the medium [21] [22]. More- 
over, there has been a lot of activity in recent years on 
implementing these dynamical mechanisms into Monte 
Carlo models of the final state parton shower [23H29] , 
including the role of decoherence and collisional broad- 
ening. Finally, the importance of decoherence of the 
partonic shower, and of collisional mechanisms has also 
been emphasized independently in [30] [31]. In the present 
rapid communication, we have combined these previously 
formulated key ideas in the dynamical picture of the 
medium as a frequency collimator of jet parton show- 
ers. We have then provided semi-quantitative estimates 
to support the view that frequency collimation can ac- 
count for the main features seen in the medium-induced 
dijet asymmetry measurements by ATLAS. We expect 
that this observation can provide some guidance for the 
further development of jet quenching models and Monte 
Carlo tools. 
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